Next  Abstract 


Lincoln  Laboratory  ASAP-2001  Workshop 


Passive  Differential  Matched-fieid  Depth 
Estimation  of  Moving  Acoustic  Sources 


Shawn  Kraut  and  Jeffrey  Krolik 


Duke  University 

Department  of  Electricai  and  Computer  Engineering 

Durham,  NC  27708 


Support  by  the  6.1  Program  of  the  Office  of  Naval  Research  Code  321  US 


niveisity 


^cho:)l  of  Engineering 


1 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

14  MAR  2001  N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Passive  Differential  Matched-field  Depth  Estimation  of  Moving  Acoustic 
Sources 

5a.  CONTRACT  NUMBER 

F19628-00-C-0002 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Shawn  Kraut;  Jeffrey  Krolik 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Duke  University  Department  of  Electrical  and  Computer  Engineering 
Durham,  NC  27708 

8.  PERFORMING  ORGANIZATION  REPORT 

NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  ADM001263  for  entire  Adaptive  Sensor  Array  Processing  Workshop., 
color  images. 

The  original  document  contains 

14.  ABSTRACT 

See  briefing  charts. 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION 

HF  APl^TRAPT 

1 8.  NUMBER  19a.  NAME  OF 

rtF  PAFIFQ  PFQP/AMQTPT  F  PFPQ/AM 

v_/Jr  /\13  o  1  1 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  jj-jj 

unclassified  unclassified  unclassified 

17 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Passive  Moving  Target  Depth  Estimation  (MTDE) 


OBJECTIVE:  To  discriminate  submerged  versus  surface  targets  by  exploiting  changes  in 
the  spatial  wavefront  at  the  array  due  to  multipath  propagation  from  a  moving  source. 

BACKGROUND: 

•  Conventional  matched-field  processors  use  a  computational  model  to  predict  the  relative 
phase  and  amplitude  between  multipath  arrivals  from  a  distant  stationary  source  and 
thus  are  very  sensitive  to  horizontal  wavenumber  differences  multiplied  by  range. 

•  Target  motion  in  classical  MFP  techniques  is  problematic  since  it  tends  to  decorrelate 
multipath  components  over  the  observation  times  used  with  stationary  source  models. 

•  Previous  work  using  moving  sources  has  attempted  to  mitigate  source  motion  effects  by 
pre-processing  so  as  to  effectively  remove  target  dynamics. 

•  Proposed  work  aims  to  exploit  target  dynamics  to  estimate  source  depth  and  range-rate 
without  requiring  the  accurate  environmental  models  required  for  range  estimation. 

•  Joint  depth-range-rate  estimation  should  achieve  robustness  to  environmental  mismatch 
since  it  depends  only  on  horizontal  wavenumber  differences  multiplied  by  the  change  in 
target  range. 
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Conventional  MFP  with  a  Vertical  or  Horizontal  Array 


•  A  snapshot  of  tilted  vertical  linear  array  (TVLA)  data  can  be  modeled  as: 

V  =  5  t/(0  )a  +r| 

where  [C/Oj)]™,  =(|)/(z^ 

0s,r^,z^  are  source  bearing,  range,  depth,  and  y  is  array  tilt. 

•  Full  3-D  range-depth-bearing  adaptive  MFP  requires  accurate  prediction  of  {k^  —  kj)r^ 
which  is  difficult  for  large  range,  sufficient  observation  time  over  which  the  source  can  be 
considered  stationary,  and  a  search  over  3  variables. 
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Some  Previous  Depth  Estimation  Approaches 


•  Averaging  a  3-D  MV  surface  over  range  may  be  computationally  intensive.  Further, 
matrix  inversion  prior  to  averaging  can  be  statistically  unstable. 

•  Matching  the  normal  mode  power  distribution  versus  hypothesized  target  depth  requires 
near  orthogonality  between  modes  at  the  array. 

•  The  MV  adaptive  beamformer  with  extended  range  constraints  (MV-ERC)  consists  of 
widening  the  range  mainlobe  so  that  bearing-depth  estimation  can  be  performed  in 
coarse  range  bands. 

•  Desensitizing  the  adaptive  beamformer  to  target  range  variation  permits  the  use  of 
longer  observation  times  for  more  stable  CSDM  estimation. 

•  The  ambiguity  surface  for  the  MV-ERC  beamformer  is  given  by: 

where  =  [1,0,...,0]''’,  where  H(r,z)  are  the  dominant  eigenvectors  of 
1  ^ 

—  Sd(r-l-  Ay^,z)d(r  -I- ,z)^  and  A^,A:  =  1,...,V defines  a  coarse  range  band  around  r. 
^k=\ 
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MV-ERC  Matched-field  Beamforming  Results 


•  Typical  simulated  ambiguity 
surface  for  8 -tonal  SWELLEX-96 
TVLA  scenario  with  SR=5800  m, 
SD=49  m.  and  SB  =  asin(0.34). 


•  Typical  real  ambiguity  surface 
for  8 -tonal  SWELLEX-96  TVLA 
event  S5  5/10/96  2332  Z.  Obs. 


Time  =  54  s,  SR=6100  m,  SD=56 
m,  and  SB=asin(0.33). 
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Fundamental  Depth  Estimation  Considerations 


Robust  bearing-depth  discrimination  without  range  could  in  principle  be  obtained  by  ^ 
treating  modal  phase  terms  as  nuisance  parameters,  i.e.  [a{rg,Zs)\i  =  ^ ^ 

Cramer-Rao  Lower  Bound  (CRLB)  on  source  depth  (left)  for  the  known  (“exact”)  versus 
unknown  (“modified”)  modal  phase  suggests  depth  estimation  without  range  possible. 

Sampled  modal  eigenfunctions  (right)  sampled  at  two  depths  illustrate  ambiguity  in 
estimating  modal  phases  jointly  with  source  depth  without  target  motion. 


£0 

SNR  IdB] 


Mode  # 
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A  Dynamical  Model  for  Passive  Depth  Estimation 


•  Idea  is  to  exploit  modal  phase  trajectory  under  a  constant  range-rate  hypothesis  in  order  to 
jointly  estimate  target  range-rate  and  depth. 

•  Letting  the  complex  range-dependent  modal  amplitudes  of  a  source  for  snapshot  k  be 
denoted  the  relative  changes  in  modal  phase  from  snapshot-to-snapshot  impose  a 
Markov  state  update: 


where  AAr)  =  diag{e^^’^  “'* '*)  and  the  additive  process  noise 

approximately  accounts  for  horizontal  wavenumber  uncertainties. 

•  The  spatial  wavefront  at  the  array,  y^,  at  narrowband  snapshot  k,  is  then  obtained  by 
taking  the  sum  of  the  normal  modes  multiplied  by  an  i.i.d.  zero-mean  Gaussian  random 
scalar, 

yk  +% 

where  <{>(2  =  diagi^iiz  ))  and  represents  additive  noise. 
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A  Recursive  Resampled  Bayesian  Estimate  for  Depth 


•  The  non-linear  depth-range-rate  estimation  problem  can  be  solved  by  representing  the  posterior  density 
function  of  the  state  by  a  set  of  random  samples,  rather  than  a  continuous  function  over  some  high 
dimensional  state  space. 

•  For  example,  suppose  at  step  k,  random  samples,  i  =  are  available  from 

p{^k-\  I  Tf •  •  • » fk-X )  •  samples,  xj^  (i)  from  p{xj^  I  yj , . ,  yj^_i )  can  be  obtained  using  these 

samples  as  input  to  the  state  equation  together  with  samples,  8  „  ,  drawn  from  its  known  distribution. 

Ok 

* 

•  The  updated  posterior  density  can  then  be  approximated  at  each  sample,  Xj^  (/),  by  forming: 

„  ^  p(yk^xpi)) 

i=i 

•  Samples,  (/),  i  =  can  now  be  obtained  by  bootstrap  resampling  times  from  the  discrete 

* 

distribution  defined  such  that  for  any  j,  (j)  =  (/)}  =  qi.  The  conditional  mean  of  the  depth 

parameter,  ,  can  then  be  estimated  by  averaging  these  bootstrap  samples. 

•  These  steps  are  repeated  for  each  range  step  to  obtain  a  recursive  estimate. 

•  For  passive  sonar,  i  \  is  zero-mean  Gaussian  with  covariance =<3^U<t>(z  )x  x  +<3f 

,  PXypPk’  k  s  s  k  k  s 

as  m  conventional  models. 
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Sequential  Importance  Sampling  Illustration 


•  Illustration  clockwise  from  upper  left  of  random  samples  from  prior,  weighting  by 
likelihood  function,  Monte  Carlo  re-sampling  from  updated  posterior,  prediction  using 
random  samples  of  state  noise. 

Sam  pi#  from  Prior,  J  1  ■  Weight  By  ukehhood, 
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Recursive  Bayesian  Passive  MTDE  Summary 


State  vector  includes  unit-magnitude  modal  coefficients  and  range-rate  with  uniform  prior.  \ 


•  State  transition  density  assumes  multiplicative  modal  phase  noise. 

•  Conditional  density  of  data  snapshot  given  state  vector  is  zero-mean  complex  Gaussian. 

•  Estimates  of  depth-range-rate  likelihood  achieved  using  sequential  importance  sampling. 


Mismatched  Range-Independent  Simulation 


•  To  illustrate  passive  moving  target  depth  estimation,  a  simple  simulation  was  performed 
using  a  normal  mode  model  of  a  15-mode  Pekeris  waveguide  with  a  100  m.  waveguide. 

•  Simulation  of  0  dB  SNR  targets,  at  2  m.  and  20  m.  depths  with  2  m/s  range-rate, 
received  at  a  water-column  spanning  23  sensor  vertical  array  using  ~50  narrowband 
snapshots  was  considered. 

•  Large  environmental  uncertainty  was  simulated  by  adding  independent  uniform  random 

variables  to  the  Pekeris  vertical  wavenumbers,  -i-M^  ,  where  =  U (-0.457r//i,0.457r//i) 

used  to  compute  both  horizontal  wavenumbers  and  modal  depth  eigenfunctions. 


Vertical  wavenumbere  lor  eigenmcdes 
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Depth  (m) 


Mismatched  Conventional  MFP  for  a  Submerged  Source 


•  Example  conventional  (aka  “Bartlett”)  ambiguity  surface  (left)  for  moving  source  at  20 
m.  depth  illustrates  extreme  ambiguity  problem. 

•  Log-histogram  (right)  of  MFP  estimates  using  100  Monte  Carlo  trials  illustrates  poor 

range  estimation  and  mediocre  depth  estimation  performance  over  =  t/(- 0.4571 //z,0.457r//z) 


BarHett  Ambiguity  Surface,  #snapshots=50 


Log-Histogram  of  Bartlett  Estimates,  target  at  20m,  #trlals=1 00 


Range  (m) 

-16.5  -16  -15.5  -15  -14.5  -14 
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Depth  (m) 


Mismatched  Passive  MTDE  for  a  Submerged  Source 


•  Example  depth-range-rate  log-likelihood  surface  (left)  for  moving  source  at  20  m.  depth 
illustrates  depth  ambiguity  with  excellent  range- rate  estimation. 

•  Log-histogram  of  MTDE  (right)  over  100  Monte  Carlo  trials  illustrating  joint  depth- 
range-rate  estimation  performance  over  =  U(-  0.457i //i,0.457r /h) 


Range  rate  (m/s) 
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Mismatched  Conventional  MFP  for  Near-Surface  Source 


•  Example  conventional  (aka  “Bartlett”)  ambiguity  surface  (left)  for  moving  source  at  2 
m.  depth  and  2  m/s  range-rate  illustrates  extreme  range  ambiguity  problem. 


•  Log-histogram  (right)  of  MFP  estimates  using  100  Monte  Carlo  trials  illustrates  poor 
range  estimation  but  good  depth  estimation  performance  over  =  U (-0.457c//z,0.457c//z) 


•  In  practice,  detection  of  a  moving  source  from  uncorrelated  surface-based  noise 
seriously  limits  the  use  of  conventional  MFP  for  depth  classification. 


Ba  fllfitt  Ambiguilv  So  dace 


Log-HislogffliTi  ot  Bartlett  Estimates,  target  at£m,  #trials=100 
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Mismatched  Passive  MTDE  for  a  Near-Surface  Source 


•  Example  depth-range-rate  log-likelihood  surface  (left)  for  moving  source  at  2  m.  depth 
illustrates  depth  estimation  comparable  to  MFP  with  excellent  range-rate  estimation. 

•  Log-histogram  of  MTDE  (right)  over  100  Monte  Carlo  trials  illustrates  good  range-rate 

estimation  and  depth  estimation  performance  over  =  U (0,0. 97r /h) 

•  Ability  of  passive  MTDE  to  discriminate  constant  range-rate  sources  may  facilitate 
detection  of  targets  at  depth  from  surface  shipping  with  different  range- rates. 
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Comparison  of  Depth  Estimation  Performance 


•  Histograms  of  MTDE  (upper  left)  vs.  conventional  (lower  left)  for  submerged  source 

with  =  U{-  0.4571  / /z, 0.4571  jh)  mismatch  indicates  moderately  improved  performance. 

•  Probability  of  correct  depth  localization  (notwithstanding  range,  range-rate,  or  predicted 
depth  ambiguity)  compares  performance  for  surface  versus  submerged  source. 


Joint  PCL  for  MFP  range-depth  estimate  versus  MTDE  range-rate -depth  estimate 
expected  to  show  significant  improvement  for  latter  approach  in  mismatched  channels. 


Histogram  of  Depth  Estimates,  MFDE,  target  at  20m,  #tnals=1 00 


PCL  (within  2  m,  amb.  inch 
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Conclusions  and  Future  Work 


•  Moving  target  depth  estimation  (MTDE)  shows  potential  as  a  classification  tool  for 
passive  discrimination  of  sources  in  the  water  column  versus  surface  ships. 

•  By  exploiting  target  motion,  MTDE  jointly  estimates  depth  and  range-rate  avoiding 
severe  range  ambiguity  problems  of  conventional  MFP  in  mismatched  conditions. 

•  Current  sequential  importance  sampling  approach  for  solving  MTDE  requires  further 
development  for  operation  at  lower  signal-to-noise  ratios. 

•  MTDE  framework  may  provide  an  alternative  detection  strategy  by  considering  bearing- 
range-rate-depth  likelihood  surfaces. 

•  Current  work  includes  evaluating  MTDE  with  real  S  ACLANT  and  SWELLEX  data. 

•  Straightforward  broadband  passive  implementation  can  be  achieved  by  incoherently 
summing  a  posteriori  probabilities  across  the  frequency  band. 
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